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We study the photon creation inside a perfectly conducting, spherical oscillating cavity. The elec- 
tromagnetic field inside the cavity is described by means of two scalar fields which satisfy Dirichlet 
and (generalized) Neumann boundary conditions. As a preliminary step, we analyze the dynamical 
Casimir effect for both scalar fields. We then consider the full electromagnetic case. The conser- 
vation of angular momentum of the electromagnetic field is also discussed, showing that photons 
inside the cavity are created in singlet states. 
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I. INTRODUCTION 



The dynamical Casimir effect consists in the generation of photons from the vacuum state of the elec- 
tromagnetic field in the presence of time-dependent boundaries or time dependent media 0, 0] ■ From the 
theoretical point of view, it is widely accepted that the most favorable scenario to observe the phenomenon 
involves a periodic time dependence, to enhance photon creation by parametric resonance. Although no 
concrete experiment has yet been performed to confirm this non-stationary Casimir effect, an experimen- 
tal verification is not out of reach, and there are several interesting proposals and ongoing experiments to 
observe it Q . 

Since the pioneering work by Moore Q , the dynamical Casimir effect for time dependent boundaries has 
been studied for different fields and geometries: scalar fields in one dimensional cavities , and in three 
dimensional rectangular 0,0, and spherical Q cavities. For the electromagnetic fieldjit was analyzed in 
three dimensional rectangular [Tol | and cylindrical cavities with time dependent length fllT | and radius |l2j| . 
The spherically symmetric situation has been considered to study quantum radiation from a time dependent 
interface between two dielectric media jl^ ■ 

In this paper we consider the quantum electromagnetic field inside a perfectly conducting, spherical cavity 
with a time dependent radius. Although a spherical conducting shell may not be appealing from an experi- 
mental point of view, it presents many interesting theoretical aspects. On the one hand, there is no classical 
electromagnetic radiation for a pulsating sphere, not even if charged. Thus it is interesting to check whether 
the quantum effect exists or not. Moreover, the angular momentum conservation implies that, if the effect 
exists, photons should be created in singlet states. On the other hand, it is also of interest to compare the 
rate of TE and TM photon creation. 

The paper is organized as follows. In Section II we describe the classical electromagnetic field inside a 
spherical cavity with time dependent radius. The TE (TM) modes are described by a scalar field satisfying 
Dirichlet (generalized Neumann) boundary conditions. For this reason, it is of interest to analyze the case of 
quantum massless scalar fields satisfying both boundary conditions, which we do in Section III. We study in 
detail the resonant case in which the cavity oscillates at twice the frequency of some field mode. In Section 
IV we quantize the electromagnetic field inside the cavity, and compute the number of TE and TM created 
photons. In Section V we discuss the conservation of the angular momentum of the electromagnetic field 
inside the cavity. Section VI contains our main conclusions. We use natural units h = c = 1. 
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II. CLASSICAL ELECTROMAGNETIC FIELD INSIDE A SPHERICAL CAVITY 

We consider a cavity bounded by a perfectly reflecting spherical shell. We will assume that the shell is at 
rest for t < 0, and that it moves following a given trajectory a(i), for < t < i/. The trajectory is prescribed 
for the problem and works as a time-dependent boundary condition for the field. Moreover, we will assume 
a non relativistic motion of the shell with a(t) = ao(l + e/(i)) with e 1 and f(t) a smooth function that 
vanishes for t < and t > tf. 

The electromagnetic field inside the cavity can be described in terms of the four vector potential = 
(ip, A). In the Coulomb gauge, V ■ A = 0, the scalar potential (p vanishes and the vector potential satisfies 
the wave equation DA = 0. 

In order to obtain the vector potential A we consider a function cj> satisfying the scalar equation □</> = 0. 
A solution of the wave equation for a given orientation of the axes must be solution for any other orientation, 
so we can obtain A starting from cf> by means of A = iLcfi = rx V<fi. This is not the more general solution for 
the electromagnetic field inside the cavity. For example, the electric field obtained from A has no component 
along r (it is a TE mode). A linearly independent solution can be obtained, however, by interchanging the 
roles of E and B. Therefore, the electromagnetic field can be written in terms of two vector potentials A TE 
and A™. Both of the them are obtained by the application of iL on solutions of the scalar wave equation, 
(j) TE and (j>™ , and satisfy the Coulomb gauge. The complete fields arc 

E = E TE + E™ = -d t A TE + V x A™ (1) 



B = B TE + B™ = V x A TE + d t A™ (2) 

The boundary conditions for the electromagnetic field on a moving interface between two media are |l4j 

(D n - Dj).n = a 
(Bjj - B/).n = 
{n x (Hj/ - Hj) + (v.n)(D„ - Dj)}.t = K.t 

{n x (E JZ - Ej) - (v.ft)(Bjj - Bj)}.t = (3) 

where n denotes the normal to the interface going from medium I to medium II, t is any unit vector 
tangential to the surface, a is the surface charge density and K the surface current. These conditions can 
be derived using the Maxwell equations in the laboratory frame or, alternatively, by performing a Lorentz 
transformation in which a given part of the interface is instantaneously at rest, and imposing there the static 
boundary conditions [l^ . 

We assume the spherical shell to be a perfect conductor, so the fields vanish in region II and the boundary 
conditions become 

B.f = {E x f + d(t)B}.t = (4) 

For the static cavity, we have a(t) = and the boundary conditions are the usual ones. In terms of the 
scalar fields these conditions read 

4> TE \ r = ao =0 [d r (r<f M )}\ r = ao = (5) 
When the shell begins to move Eq. J3J implies 

4> TE \ r =a(t) = {(d r + a(t)d t )r(f>™ }\ r=a(t) = (6) 

From this discussion we see that the behavior of the vector potential A TE (A™)is related to the problem 
of a scalar field subjected to Dirichlet (generalized Neumann) boundary conditions. The description in terms 
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of independent TE and TM fields is possible due to the particular geometry we are considering. Indeed, using 
the above definitions and boundary conditions it is easy to check that no mixed terms appear in Maxwell's 
Hamiltonian. 

The two scalar functions 4> TE and (f)™ arc known as Debye potentials 0|. These functions are related 
with the Hertz potentials IT e and II m used in Ref. by 



n e = ^r n m = ^r (7) 

For the sake of simplicity, instead of dealing with the full electromagnetic case, we will first study the problem 
of scalar fields. 



III. QUANTUM SCALAR FIELDS INSIDE A SPHERICAL SHELL 

Let us consider a masslcss scalar field <fi(r, t) inside a spherical cavity described by the Lagrangian C 

c = \d^d^ (8) 

The field operator can be written in terms of creation and annihilation operators as 

< M r ' t) = tt MmWm(r, t) + h.c. (9) 
k£m 

where the mode functions <PkZm{ v i t) form a complete orthonormal set of solutions of the wave equation. 

A. Static cavity 

When KOwe have a static spherical shell of radius a and the field modes are given by 

fk£rn(r,t) = (j} k l m {r)—== (10) 

where 

4>k£ m {r) =Cke m 3t(u£kr)Ye m (9,4>) (11) 

are the eigenfunctions of the Laplacian with eigenvalues — a->f fc , ji are the spherical Besscl functions and Yi m 
the spherical harmonics. The normalization constants Ckim depend on the boundary conditions satisfied by 
the field. 

We will consider two scalar fields, one satisfying Dirichlct boundary conditions and the other one satisfying 
generalized Neumann boundary conditions. It is worth to note that the scalar fields considered here are not 
exactly the Debye potentials of the electromagnetic field (the Lagrangian in Eq.JSJ is not the Maxwell 
Lagrangian expressed in terms of the Debye potentials). In spite of this, we will denote them by (f) TE and 
(ft™ to stress the boundary conditions that they satisfy. 

For the scalar field (j) TE the normalization constant is given by 



1 

\]'e{jtk) 



Ckim = \l -3 -rr s (12) 



with jik the k-th zero for the spherical Bessel function ji(x). The frequencies of the modes are uiik = ^ 



4 



On the other hand, for the scalar field <j>™ we have 

Cklm = {^W^)7^W+% (13) 

where is the k-th zero of {d x [xje(x)]} = 0. The TM frequencies are given by ugh = 

The operators a^ m and a^" m create and annihilate particles with well defined energy, total and z- 
component of the angular momentum £ and m respectively. They correspond to the particle notion in 
the "in" region (t < 0). 



B. Moving cavity 

1. Dirichlet Boundary Condition 

When the radius of the cavity depends on time, the modes of <p TE can be expanded in terms of an 
instantaneous basis 

t) = E Q P %E<t>pL(r, «(*)) (14) 
p 

with 



^ E m (r,a(t)) = JirtS)Ii»(M) (15) 

V a ( tJftOtp) a w 

Because of the spherical symmetry, in the expansion of the mode (r, t) it is enough to use the functions 
<fipfo n with the same values of £ and m, i.e. we only mix the first quantum number. Although the coefficients 

Qp%E depend on the angular momentum £ (see Eas. (|16ll7H below), in order to keep the notation as simple 
as possible we do not write this explicitly. 

The initial conditions for the coefficients Qp TE (t) are 



Ql%(t = 0) = -^== = 0) = -y ^4, 



These conditions ensure that, as long as a(t) and a(t) are continuous at t = 0, each field mode and its time 
derivative are also continuous functions. 

The expansion in Eq. (|14|l for the field modes must be a solution of the wave equation. Taking into account 
that at each time the functions </>J^(r, a(i)) form a complete and orthonormal set, the wave equation is 

equivalent to the following set of coupled equations for Q„ TE (t) 

QnjrEii) + Mt)?Q ( n %(t) = -2X(t)J2Qpi E (Hn 

p 

~Ht)Y,Qp%(Hn+0^) (17) 



where \(t) = and 



r 2rr r n r a(t) f)fh TE 
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The coefficients gi n can be computed explicitly using that, for k ^ k', 



r 2 Je (kr)3e(k>r)dr = _ fc2 {k M (k' a^ka) - k' M {ka)j',{k' a)} (19) 



The result is 



i i f if p = n 

9 P n = -9n P = \ 2 ff"-»p if p ^ n . (20) 

2. Neumann Boundary Condition 

We now consider a scalar field satisfying the generalized Neumann boundary condition on the surface of 
the shell. To satisfy the boundary conditions for t > we will expand the mode functions with respect to 
an instantaneous basis, as we did for the case of Dirichlet boundary condition. However, in this case the 
choice for the instantaneous basis is not so easy as before, because the boundary condition on the moving 
shell given in Eq.© involves a time derivative of the field. 

The instantaneous basis can be obtained by means of a change of variables in the (t, r) plane |1C( . provided 
that in the new variables (77, £) the boundary condition is the usual one (i.e. no time derivative of the field) 

{dd&™(V,Z,0,<P)}}k=l(r l ) = (21) 

where £(77) is the value of the coordinate £ on the moving spherical mirror. 

We define the line ;/ = const to be a slight modification of the line t = const, in such a way that it is 
orthogonal to the worldline of the mirror (t, a(t)r) at r = a(t). The coordinate £ is defined as the distance 
from r — to r on the line 77 = const. Explicitly 

V = t + g(r,t) 



£ = f r dr\l + ( 22 ) 

where g(r, t) = O(e) and therefore £ = r + 0(e 2 ) and l(n) = a(t) + 0(e 2 ). 

In order to ensure the orthogonality between the line 77 = const and the world line of the mirror, we impose 

m 

g(a(t), t) = d r g(r, t)\ r=a[t) = -a(t) (23) 

The function g(r, t) is of course not unique. It can be expressed as g(r, t) = a(t)a(t)v(r / a(t)) , where v(l) = 
and i>'(l) = — 1 (the prime denotes derivation with respect to the argument). There are many solutions to 
these conditions, implying a freedom for selecting the instantaneous basis. However, physical quantities like 
the number of created particles, the energy density inside the cavity and the angular momentum of the field 
are independent of the particular choice of g(r,t) [Tol |. 

In the new coordinates, the instantaneous basis is the trivial one 



fklm 



Returning to the (t, r) variables, each field mode can be expanded as follows 

^™{t, r, 9, <p) = J2[Qp^TM(t) + Q { p %(t)9(r, tM™(r, a(t)) + 0(e 2 ) (25) 
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with 



a 3 (t) JiM ^/ K 2 £p -£(£+!) 



As in the Dirichlet case, the coefficients Q^xm depend on the number £, but we do not write the dependence 
explicitly. Assuming that ait) and d(t) are continuous at t = 0, and that the initial acceleration satisfies 
d(0) = 0(e 2 ), the initial conditions for Q p tm$) arc * nc same as those for Qp% E (t), Ea. pti|) . The equation 
of motion for \, M (t) is 



Qi%i(t) + Nn(t)] 2 Ql%(«) = -2A(t) E P Q^Mflpn - A W E P Q^wi 



-2a 2 (t)A(t)E p Q^M^ - E p Q^M[4n A (*) fl2 (*) - WHn 

)(*0 ,< 
»p,TM pn 

where the coefficients Sp„ , and c/p„ are given by 

/•2-7T />7T ra(t) 

'o Jo Jo 



c — 
pn 



-A(t)a 2 (i) V d*QWs* + 0{e 2 ) (27) 



^ OT ^(r,a(i))C(r,«(t)) (28) 



»£. = JcT So Io a(t) d 3 xa 2 (t){[d 2 r r v - (L, tp )*v]<t%&(T, a(t))tf™* (r, a(t)) 

+^ r t;fl r [r^(r >a (t))]^(r >a (t))} (29) 



if p =n 

fpn = S " / K 2 (30) 



C. Creation of particles 

We are interested in the number of particles created inside the cavity, so it is natural to look for harmonic 
oscillations of the shell that could enhance that number by means of resonance effects for some specific 
external frequencies f2. So we study the following trajectory 

a(t) =a (l + esin(fM)) (31) 

Let us first consider the Dirichlet scalar field (j> TE . When t > tf, ("out" region), the radius returns to its 
initial value ag, the right hand side in Ea. ll7|) vanishes and the solution is 



QuJte^ > */) — A tn>TE e wtnt + B inTE e (32) 

constant coefficients to be deb 
In these coefficients we write the dependence on £ explicitly. 



where A^ TE and B^ TE are constant coefficients to be determined by the continuity conditions at t = tf. 
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For t > tf we can define a new set of operators a^m anc ^ a ktmi associated with the particle notion in the 
"out" region . The "in" and "out" operators are connected by means of the Bogoliubov transformation 



^„TE a kim + ( - 1 ) ™ A* TE a kt-m\ ^ 2uJ Hn (33) 



The number of "out" particles in the mode (n,£,m) is given by 



< ATrUm >=< l n\a° n fla° n f m \0 m >= 2uj, n |4StbI 2 (34) 

k 

To obtain the coefficients A$ TE and Bgj TE we must solve Ea. (|17|l for the coefficients Q„r#(i)- They are 
of the same form as those that describe the modes of a scalar field satisfying Dirichlct boundary conditions 
in a three-dimensional rectangular cavity Q, and can be solved using Multiple Scale Analysis (MSA), see 
for example, Refs.0,0]. We will review here the main results and include the details in the Appendix. 

The solution for the coefficients A^ TE and B^ TE depends on the relation between the external frequency 
n and the natural frequencies of the field in the cavity. There is parametric resonance when the external 
frequency CI equals the sum of two cigenfrcquencics of the cavity with the same angular momentum Q = 
tU£ n + LOi p . For simplicity, in what follows we will analyze the particular case f2 = 2uj£ n . Using MSA, one can 
show that the modes (£, n) and (£, q) are coupled if any of the following conditions is satisfied 

^ = Win — ^iq f& = ~ ^fn + ^>tq (35) 

For $ TB the frequencies are given by the zeros of the spherical Bessel functions (u>g n = and so, the 
spectrum is qualitatively different depending on the value of the number I. For I = the spectrum is 
equidistant. In the particular case = 2cjon, the set of equations for the coefficients Aq^ te and B^ TE 
corresponds to that of a one-dimensional cavity excited with twice the lowest eigenfrequency. The equations 
are coupled, the number of particles grows linearly in time, and the energy increases exponentially jg. 

On the other hand, for t ^ the spectrum is not equidistant [l^, and one can check that there is no 
coupling between modes. Ea. (|17|) reduces to the Mathieu equation for the modes with frequency u>£ n . The 
Bogoliubov coefficients A^ TE and B^ TE and the number of created particles grows exponentially. 

For the case of the scalar field <j>™ the situation is similar. The set of Eas. l|27|) has the same form 
as the set that describes the modes of a scalar field in a three-dimensional rectangular cavity satisfying 
generalized Neumann boundary conditions ^(j; an d can be solved again using MSA. As before, the spectrum 
is equidistant for I — and non-equidistant for i ^ 0. 

In the resonant case = 2ojln (L ^ 0), for both TE and TM modes the number of created particles is 
given by (see Appendix) 

< AW >=< in \a™2 m a%i m \0 in >= sinh 2 ( 7 et f ) (36) 

and 

< Mnl >= E < Ma^l^l JOm >= (2i + 1) Binh 2 (7ct / ) (37) 

m 

The constant 7 determines the rate of growth. For TE modes we have 

1 TE = (38) 
2a 



while for TM modes 

1 



L(L+l) 



1 -2a Ql -^tll [6J) 



LN 



The case I = is qualitatively different and, as stated above, equivalent to the one dimensional dynamical 
Casimir effect. However, as we will see in the next section, these modes are absent for the electromagnetic 
field. 
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IV. THE ELECTROMAGNETIC FIELD 

In Section II we showed that the electromagnetic field inside the spherical cavity can be described in terms 
of two vector potentials. They can be expanded in terms of creation and annihilation operators as 

A TE (r , t) = < m A-lL (*,*) + h.c. (40) 

kirn 

A™(r, t)=J2 4" m A™ ( r , t) + h.c. (41) 

kirn 

The modes A^ E n (r,t) and A^^(r, t) can be obtained from the modes of two scalars fields through the 
application of the operator iL = rx V . This operator acts only on the angular part of the scalar modes, so 
we introduce the vectorial spherical harmonics ~K.( rn (9, <j)) 

Xt m (9, <f) = —±== Y em (0, <f>) (42) 



(the additional factor y/£(£ + 1) is needed for normalization) . Therefore the modes for the vector potentials 
are given by 



7 e Utkj- 

A££(r,t < 0) = vf^T ^Vi) ^^^^^ 



for the static cavity and 



A2£.(r,t > o) = E P y^^Q^M^xU^) 

A^(r,t > 0) = E P V^^7^=Wf) [0 ^ W + 5M)Q^ M «b^)XUM) (44) 

for the moving cavity. It is worth to remark that, as LY"oo = 0, there is no monopolar term in the expansions 
for the electromagnetic field. 

The dynamical evolution of the TE (TM) modes is that of the modes of 4> TE {4>™) with I ^ 0. As a 
consequence, the number of created photons in each mode equals the number of created particles of the 
corresponding scalar field. If we consider again the parametric resonant case fi = 2loln, the number of 
photons grows exponentially and is given by the Eas. (|36l38|) in the case of TE modes and by Eqs. (|36l39fl in 
the case of TM modes. It is interesting to note that it is in general not possible to excite at the same time 
both a TE and a TM mode (for that to be possible one should have loJ^ + uifjf = loJ,^ 1 + , which is 
not satisfied). 

V. ANGULAR MOMENTUM 

In this section we discuss the conservation of the angular momentum of the electromagnetic field inside 
the spherical cavity, showing that photons are created in singlet states. As we start in the vacuum state 
of the electromagnetic field, the average value of the angular momentum is initially zero. The oscillations 
of the radius of the cavity does not break the symmetry under rotations, so the angular momentum must 
vanish at all times. 
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The angular momentum of the electromagnetic field is 

L = f r x (E x B)d 3 x = L TB + L™ (45) 

J sphere 

As in resonant situations it is possible to produce either TE or TM photons, in what follows we will consider 

rp rp rp 71 f 

only one of the two polarizations, without specifying which one. Both L and L are of the form 

L = (LuIjjjLz) = J2ke7n(^i^tm-l a klm a k£m-l + ^trn a ktm- 1 a kim\ J 

=l\r + n m ^ n in - C~ n in>{ n in 1 m n m ^ n in \ I A.R\ 

2 ' lm a ke.m+l a klm ^ (.m a kf,m-l a klm\^ " La ke.m a klm) \ w I 

with 

C+ n = y/y-mW + m + l) Cj m = y/(e + m)(e-m+l) (47) 

The expression of L in terms of the " in" creation and annihilation operators is valid for all times. Further- 
more, as we work in Heisenberg picture, the state of the electromagnetic field is always the "in." vacuum, so 
we have 

< L >=< m |L|0 m >= Vf (48) 



However, in the "out" region, when the radius of the cavity returns to its original value and remains at 
rest, we could expand the field in terms of the new set of the "out" creation and annihilation operators. 
Therefore the angular momentum of the field can be written as in Ea. (|46J) but changing the "in" operators 
by the "out" ones. As the number of "out" particles is different from zero, photons must be created forming 



singlet states, in such a way that the angular momentum of the field remains null. 
The state of the field |0j„ > can be written as a linear combination of "out" states 

|0 m >= a\0 O ut > + X] a nlraa°^\Qout > + X/ X/ a ntm,n'l'm' a°"„ Vf'm' \®»Ut > + (49) 

n£m nlm n' £' m' 

In the parametric case 57 = 2cj£/v, the relation between the "in" and "out" operators is, for I = L and 

n = N, 

a° N u l m = [4'4 m + (-l) m (4$r<#L JV2^ (50) 

< L l = [(4$)*<#L + (-ir4>'^- (si) 

The "in" and "out" operators coincide a°™^ = a^™ m when £ ^ L or n ^ N. Therefore, we rewrite Ea. (T4*9l) 
in the following way 

|0j„ >= Q(|0 o „ t > +y^Qim|loMt >ro +^^am,m'|loMt >m |lo«t >m> + (52) 

m m m' 

where we omitted the subindexes TV and L. 



In order to find the a coefficients we apply a destruction "in" operator on both sides of Eq. (|52|l . The 
left hand side gives zero. On the right hand side we write the "in" operator in terms of the "out" operators, 
by inverting Eas. H5U[l - i|51[l . Doing that we obtain a linear combination of orthogonal states which equals 
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zero, so the coefficient of each state must vanish. In this way we get the equations that determine the a 
coefficients. In the particular case of ft = 2usln with L=l, the state of the field |0i„ > can be written as 

|0 m >= A{\0 out > -C[\l out >i |l out >-i -^\2out >o] + 
C 2 [|2 OMt >i \2 out >_! — -i=|l out > x \2 out > 1 1^ >_i + ^^|4 out > ] — 
C 3 [|3 ou t >i \3 ut >-i — ^|2 ou t >i |2 ou t >o \2 out >_i + 

Sellout >i Kout >o lUt >-i -%§out >o] + ■■•} (53) 
where A is the normalization constant of the state and 

C = -tanh(7ei/) . (54) 

The constant 7 is given by Ea. H38|) for TE photons and by Eq. <|39|) for TM photons, both with L = 1. Each 
state between brackets in Eq. (|53|l is an cigenstate of L 2 and L z with eigenvalue equal to zero, as expected. 

VI. CONCLUSIONS 

In this paper we have computed the resonant photon creation inside a spherical oscillating cavity taking 
into account the vector nature of the electromagnetic field. We described the TE and TM modes of the 
electromagnetic field by masslcss scalar fields satisfying Dirichlet and generalized Neumann boundary con- 
ditions. We first studied the creation of particles for these scalar fields, and then showed that the number 
of created photons in each mode (TE or TM) equals the number of created particles of the corresponding 
scalar field. We used MSA to take into account resonant effects at long times, and found an exponential 
growth in the number of created photons. Previous works studied the case of a Dirichlet scalar field only in 
the short time limit eflt << 1 0. Our results are consistent with those in this limit. 

The spectrum of the scalar fields is equidistant for I — and non equidistant for £ ^ 0. When the external 
frequency is chosen to produce parametric resonance in modes with £ = 0, an infinite number of modes are 
excited. The problem becomes equivalent to the one dimensional dynamical Casimir effect. The modes with 
£ =/= are non equidistant, and there is no mode coupling. In the parametric resonance case f2 = 2a^„, 
the number of motion-induced particles grows exponentially in the particular modes corresponding to these 
values of n and £, for all possible values of m. 

The eigenfrequencies of the TE and TM modes of the electromagnetic field are the same as the eigen- 
frequencies of cj) TE and cj>™ (the only difference is that in the electromagnetic case the £ = modes are 
absent). The growth rate is different for TE and TM modes. In cubic cavities, where the frequencies of the 
TE and TM modes are equal, they satisfy 01 

1 TE < I 1™ > I (55) 
For spherical cavities, the TE and TM modes have different eigenfrequencies. The growth rates satisfy 

, .TE . TM 

1 TE = — 1™ > — (56) 

Therefore the ratio between the rate of growth and the eigenffequency is, as in cubic cavities, larger for 
TM modes than for TE modes. In other words, the functional dependence of the rate of growth with the 
cigcnfrequcncy is different for TE and TM modes. 

Finally, we considered the conservation of the angular momentum of the electromagnetic field in the 
oscillating cavity. As the initial state of the field is the vacuum, and the motion of the shell does not break 
the symmetry of rotation, the angular momentum of the field must be zero for all time. Therefore, particles 
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must be created in singlet states. For the particular case Q = 2luin we wrote explicitly the "zn" vacuum 
state as a linear combination of "out" singlet states. 

In classical electromagnctism there is no electromagnetic radiation with spherical symmetry. Therefore it 
was not completely obvious that the dynamical Casimir effect would occur in a spherical oscillating cavity. 
A mechanical analogue of the one dimensional dynamical Casimir effect is useful to illustrate this point 
[2Cj . At the classical level it is not possible to amplify transversal oscillations on a string by changing its 
length, unless there is an initial classical wave on it [20ll2ll|. At the quantum level, the initial conditions of 
the modes are non trivial due to Heisenberg uncertainty principle, and therefore it is possible to excite the 
system even starting from the ground state. Analogously, in the spherically symmetric case, the dynamical 
Casimir effect is non trivial because all modes with £ ^ have non vanishing quantum fluctuations. 
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VIII. APPENDIX 

In this Appendix we solve Eas. i|17fl for an oscillating motion of the shell given by a(t) = ao(l + esin(fM)). 
For small amplitudes (e <C 1), Eqs. I|17|) take the form 

Qi% E (t) + MQ { nJr E (t) = 2esm(nt)[Lu en } 2 Q ( n % E (t) 
+2eO cos(Ot) £ P Q ( p % E (t)g e pn - etf sin(fit) £ p Q p % E (t)g e pn 

+eW cos 2 (m) J2 pN Q^UHn^n + 0(e 2 ) (57) 

It is well known that a naive perturbativc solution of these equations in powers of e breaks down after a 
short amount of time, of order (eft)" 1 . This happens for those particular values of the external frequency 
n such that there is a resonant coupling with the eigenfrequencies of the static cavity. In this situation, to 
find a solution valid for longer times (of order (e _2 51 -1 )) we use the MSA technique 0,^3- We introduce a 
second time scale r = et and expand Qn^ TE {t) as follows 

Qi%E = Q£rk*. r) + eQ^Ut, r) + 0(e 2 ) (58) 

The initial conditions read 

Q^Ut = 0) = -J= QK(* - 0) = -tffi* (59) 

Replacing Eq.|J5SJ into Eg. (|57fl . to zeroth order in e we get the equation of an harmonic oscillator. The 
solution is 

Q^Ut, r) = Ag TB (r)e^»* + s£ TB (r)e— * (60) 
and using the initial conditions it follows that 

dtn, TE (T = 0) = B £Ite(t = 0) = -^=S kn (61) 
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To first order in e we obtain 

dttQn,T E (t,r) + k„] 2 Qi*£k*,T) = -2d Tt Q { n k f E (t,r) + 2sm(nt)[cu en } 2 Q { n k ^ E (t 7 r) 

+2Q cos(Ot) E p - tt 2 sin(fit) £ p 0^(*. (62) 

The functions A^ te {t) and B^ te (t) are obtained by imposing that no secular terms appear in the 

equation for Q„!p E (t, r), that is, any term with a time dependence of the form e± luJikt in the right-hand side 
of Ea. l|62|l must vanish. We get 

dr4 k l TE = -^4n,T E (r)S(n - 2u ln ) + E^„(§ - ^ p )B^ TE{T)g ? pnl g- 

(a* - q)5{Q + u in - ^)]ig T£ (T) S ;^ (63) 



H, - q)6(n + ut n - ^ P )]B ( ^ TE g l pn ^: (64) 

The equations above lead to resonant behavior if 

Q = 2iot n or Q. = ujin + u>£ p (65) 

These are the resonant conditions. The modes (£, n) and (£, q) are coupled if any of the following conditions 
is satisfied 

il = LUfn — L>e q fi = —Win + UJ£ q (66) 

These are the intermode coupling conditions. 

The eigenfrequencies are determined by the zeros of the spherical Bessel functions (u>g n = ^ a ). For £ = 
the spectrum is equidistant and therefore the intermode coupling conditions are satisfied. However, as for 
the electromagnetic field there are no modes with I = 0, in what follows we will consider only the case £ 7^ 0, 
in which the spectrum is not equidistant |l9j , and one can check that the intermode coupling conditions are 
not satisfied. 

For the particular case ft = 2ujln, Eas. (|63l64[) become 



d rB[% E = -^A^JeLSnN (68) 



Using the initial conditions H61[l we obtain, for £ = L and n — N 



41 rE = ^sinh(^r) 



(70) 
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and for I ^ L 

A (k) _ n R (fe) _ &kn_ / 71 s 

A tn,TE — U n £n,TE ™ ^/2uJ^ ^ ' 

If the motion of the shell ends at t = tf, the number of created particles is given by 

< M NLm >=< OinKHa^iJOin >= smh 2 (^et f ) (72) 

< M NL >=£< Vin\<L<*NL m \Qin >= (2L + 1) sintf^rf/) (73) 

m 

In this Appendix we only considered the case of the scalar field (/> TE , associated to the TE modes. The 
Eqs.(j23) for the scalar field <p™ can also be solved using MSA, following the procedure described here. 
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